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Abstract 

I show that the zeros of the chromatic polynomials Pg(q) f° r the general- 
ized theta graphs 0( s ' p ) are, taken together, dense in the whole complex plane 
with the possible exception of the disc \q — 1| < 1. The same holds for their 
dichromatic polynomials (alias Tutte polynomials, alias Potts-model partition 
functions) Zc(q,v) outside the disc \q + v\ < \v\. An immediate corollary is 
that the chromatic roots of not-necessarily-planar graphs are dense in the whole 
complex plane. The main technical tool in the proof of these results is the 
Beraha-Kahane-Weiss theorem on the limit sets of zeros for certain sequences 
of analytic functions, for which I give a new and simpler proof. 
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1 Introduction 



1.1 (Di)chromatic polynomials and the Potts model 

The polynomials studied in this paper arise independently in graph theory and 
in statistical mechanics. It is appropriate, therefore, to begin by explaining each of 
these contexts. Specialists in these fields are warned that they will find at least one 
(and perhaps both) of these summaries excruciatingly boring; they can skip them. 

Let G = (V, E) be a finite undirected graph 1 with vertex set V and edge set E. 
For each positive integer q, let Pg{<i) be the number of ways that the vertices of G can 
be assigned "colors" from the set {1,2, . . . , q} in such a way that adjacent vertices 
always receive different colors. It is not hard to show (see below) that Pg{<i) is the 
restriction to Z + of a polynomial in q. This (obviously unique) polynomial is called 
the chromatic polynomial of G, and can be taken as the definition of Pg{<i) for 
arbitrary real or complex values of q. 2 

The chromatic polynomial was introduced in 1912 by Birkhoff [IE] . The original 
hope was that study of the real or complex zeros of Pg{(i) might lead to an analytic 
proof of the Four-Color Conjecture [THEI], which states that -Pg(4) > for all loopless 
planar graphs G. To date this hope has not been realized, although combinatoric 
proofs of the Four-Color Theorem have been found [Hl2*ll31 l8*%|lll4| . Even so, the zeros 
of Pg{q) are interesting in their own right and have been extensively studied. Most of 
the available theorems concern real zeros [T71 IT27H IT331 ITM1 IT351 loTl IT?^ ITT51 1321 ITTE] . 
but there have been some theorems on complex zeros O El El 11281 12H 1221 1221 12H 
123 II 1 H I112| |2Hj as well as numerical studies and computations for special families 

[H^i fmi fT^i ei ihi ihtti nrrm [tttti [ttt^i frnm rnni frrMi n?ni frm nryHi [rnni nTT7i [ttthi tttm 
E3 H3 [[HI mi ED E21 1331 El ESI EH1 E3EHllISl[IllIiniEHlE21Eniinn]. 

A more general polynomial can be obtained as follows: Assign to each edge e E E 
a real or complex weight v e . Then define 



where the sum runs over all maps o: V — > {1, 2, . . . , q}, the 5 is the Kronecker delta, 
and x±(e) , X2(e) G V are the two endpoints of the edge e (in arbitrary order). It is not 
hard to show (see below) that Zc(q, {v e }) is the restriction to q G Z + of a polynomial 
in q and {v e }. If we take v e = —1 for all e, this reduces to the chromatic polynomial. 
If we take v e = v for all e, this defines a two- variable polynomial Zc{q,v) that was 
introduced implicitly by Whitney [1301 E3D E32J and explicitly by Tutte [UHl ITP3] : 

1 In this paper a "graph" is allowed to have loops and/or multiple edges unless explicitly stated 
otherwise. 

2 See PH103] for excellent reviews on chromatic polynomials, and for an extensive bibliography. 
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it is known variously (modulo trivial changes of variable) as the dichromatic poly- 
nomial, the dichromate, the Whitney rank function or the Tutte polynomial 

[T231ITT]. 3 

In statistical mechanics, is known as the partition function of the q-state 

Potts model. In the Potts model [HOJ [T3BJ HUH], an "atom" (or "spin") at site 
x G V can exist in any one of q different states (where q is an integer > 1). The 
energy of a configuration is the sum, over all edges e G E, of if the spins at the 
two endpoints of that edge are unequal and — J e if they are equal. The Boltzmann 
weight of a configuration is then e~@ H , where H is the energy of the configuration 
and (3 > is the inverse temperature. The partition function is the sum, over all 
configurations, of their Boltzmann weights. Clearly this is just a rephrasing of (jl.lj) . 
with v e = e /3Je — 1. A coupling J e (or v e ) is called ferromagnetic if J e > (v e > 0) 
and antiferromagnetic if — oo < J e < (— 1 < v e < 0). 

To see that Z G (q, {v e }) is indeed a polynomial in its arguments (with coefficients 
that are in fact or 1), we proceed as follows: In (jl.lj) . expand out the product over 
e G E, and let E' C E be the set of edges for which the term v e S(a(xi(e)), cr(x2(e))) 
is taken. Now perform the sum over configurations a: in each connected component 
of the subgraph (V, E') the spin value a(x) must be constant, and there are no other 
constraints. Therefore, 

Z G (q,{v e }) = J2^ E/) U V ^ (1-2) 

E'CE e£E' 

where k(E') is the number of connected components (including isolated vertices) 
in the subgraph (V,E'). The expansion (jl.2j) was discovered by Birkhoff [TSJ and 
Whitney j!3Uj for the special case v e = — 1 (see also Tutte j!18j 1119] ); in its general 
form it is due to Fortuin and Kasteleyn [H21 HE] (see also We take (jl.2j) as the 

definition of Z G (q, {v e }) for arbitrary complex q and {v e }. 

In statistical mechanics, a very important role is played by the complex zeros of 
the partition function. This arises as follows j!41j : Statistical physicists are inter- 
ested in phase transitions , namely in points where one or more physical quantities 
(e.g. the energy or the magnetization) depend nonanalytically (in many cases even 
discontinuously) on one or more control parameters (e.g. the temperature or the mag- 
netic field) . Now, such nonanalyticity is manifestly impossible in (II. lj) / (jl.2j) for any 

3 The Tutte polynomial T G (x,y) is conventionally defined as |1291 p. 45] [H] PP- 73, 101] 

T G (x,y) = {x-l) k ^- k ^{y-l)\ E '\ +k ^-\ v \ 

E'CE 

where k(E') is the number of connected components in the subgraph (V, E'). Comparison with l|1.2|) 
below yields 

T G (x,y) = (x-l)- k ^(y-l)-^Z G ((x-l)(y-l), y-l) . 
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finite graph G. Rather, phase transitions arise only in the infinite-volume limit. That 
is, we consider some countably infinite graph G^ = (Voodoo) — usually a regular 
lattice, such as Z d with nearest-neighbor edges — and an increasing sequence of finite 
subgraphs G n = (V n ,E n ). It can then be shown (under modest hypotheses on the 
G n ) that the (limiting) free energy per unit volume 

f Gac (q,v) = lim \Vn\~ 1 log Z Gn (q,v) (1.3) 

n~ >oo 

exists for all nondegenerate physical values of the parameters 4 , namely either 

(a) q integer > 1 and — 1 < v < oo [using see e.g. |SSJ Section 

1.2]] 

or (b) q real > and < v < oo [using (jl.2j) : see [53J Theorem 4.1] and 

152 EH]- 

This limit fc^iq^v) is in general a continuous function of v; but it can fail to be a 
real-analytic function of v, because complex singularities of log Zc n {q, v ) — namely, 
complex zeros of Zc„{q, v) — can approach the real axis in the limit n — > oo. There- 
fore, the possible points of physical phase transitions are precisely the real limit points 
of such complex zeros (see Theorem 13.1(1 . As a result, theorems that constrain the 
possible location of complex zeros of the partition function are of great interest. In 
particular, theorems guaranteeing that a certain complex domain is free of zeros are 
often known as Lee-Yang theorems. 5 

In summary, both graph theorists and statistical mechanicians are interested in 
the zeros of Zg(q,v). Graph theorists most often fix v = — 1 and look for zeros in 
the complex g-plane, while statistical mechanicians most often fix q real > (usually 
but not always an integer) and look for zeros in the complex f-plane. But these 
inclinations are not hard-and-fast; both groups have seen the value of investigating 
the more general two- variable problem. The analysis in this paper will, in fact, be an 
illustration of the value of doing so, even if one is ultimately interested in a specific 
one-variable specialization of Zc(q,v). 

1.2 A Lee- Yang theorem for chromatic polynomials? 

Let me now review some known facts about the real zeros of the chromatic poly- 
nomial Pc{q), iu order to motivate some conjectures concerning the complex zeros: 

4 Here "physical" means that the weights are nonnegative, so that the model has a probabilistic 
interpretation; and "nondegenerate" means that we exclude the limiting cases v = — 1 in (a) and 
q = in (b), which cause difficulties due to the existence of configurations having zero weight. 

5 The first such theorem, concerning the behavior of the ferromagnetic Ising model at complex 
magnetic field, was proven by Lee and Yang [HI] in 1952. A partial bibliography (through 1980) of 
generalizations of this result can be found in |f>5| . 
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1) It is not hard to show that for any loopless graph G with n vertices, (— l) n PG{q) > 
for real q < [HI]. It is then natural to ask whether the absence of negative real 
zeros might be the tip of the iceberg of a Lee- Yang theorem: that is, might there exist 
a complex domain D containing (— oo, 0) that is zero-free for all Pg? One's first guess 
is that the half-plane Re q < might be zero- free [33] . This turns out to be false: for 
about a decade, examples have been known of loopless graphs G that have chromatic 
roots with slightly negative real part EH1 ESJ QUI DM 121 123 E2] ; and examples 
have very recently been constructed with arbitrarily negative real part jlU5l equation 
(3.12)]. Nevertheless, it is not ruled out that some smaller domain D D (— oo, 0) 
might be zero-free. 

2) For any loopless planar graph G, Birkhoff and Lewis [17J proved in 1946 that 
Pc(q) > for real q > 5; 6 we now know that Pg(4) > [3 [EE|; and it is very likely 
true (though not yet proven as far as I know) that Pg(q) > also for 4 < q < 5. Thus, 
it is natural to conjecture that might exist a complex domain D containing (4, oo) [or 
(5, oo)] that is zero-free for all planar Pq. One's first guess might be that Reg > 4 
works. This again turns out to be false: examples are known of loopless planar graphs 
G that have chromatic roots with real part as large as ~ 4.2 EDI - Nevertheless, it 
is not ruled out that some smaller domain D D (4, oo) might be zero-free. 

As with most of my conjectures, these two are false; but what is interesting is that 
they are utterly, spectacularly false, for I can prove: 

Theorem 1.1 There is a countably infinite family of planar (in fact, series-parallel) 
graphs whose chromatic roots are, taken together, dense in the entire complex q-plane 
with the possible exception of the disc \q — 1| < 1. 

As far as I know, it was until now an open question whether the closure of the set of 
all chromatic roots (of all graphs taken together) even has nonzero two-dimensional 
Lebesgue measure. Theorem II. II answers this question in a most spectacular way. 

The graphs arising in Theorem II .11 are "generalized theta graphs" 0( s ' p ) obtained 
by parallel-connecting p chains each of which has s edges in series. 7 Theorem 11.11 is 
in fact a corollary of a more general result for the two- variable polynomials Zc{q,v): 

Theorem 1.2 Fix complex numbers qo,Vo satisfying \vq\ < \qo + Vq\. Then, for each 
e > 0, there exist s G N and a map p : N fl [s , oo) — > N such that for all s > s and 
V > Po{s): 

6 See also Theorem 1] and |115l Theorem 3.1 ff.] for alternative proofs of a more general 

result. 

7 More generally, the generalized theta graph Sl) ... lS consists of endvertices x,y connected by 
p internally disjoint paths of lengths Sx,...,s p > 1 |2(i| . The graphs arising in Theorem 11.11 thus 
correspond to the special case s\ = . . . = s p = s. See Section I2"3l below for a computation of the 
Potts-model partition function for an arbitrary Sl! ..., s . 
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(a) If Vq 7^ 0, then Zq{ 8 , p ){ • , v ) has a zero in the disc \q — q \ < e. 

(b) ZQ(s, P )(qo, ■ ) has a zero in the disc \v — Vq\ < e. 

(Setting vq = — 1, Theorem II .2f a) implies Theorem 11.11 ) 

As Jason Brown pointed out to me (why didn't I notice it myself?), an immediate 
corollary of Theorem 11.11 is: 

Corollary 1.3 There is a countably infinite family of (not-necessarily-planar) graphs 
whose chromatic roots are, taken together, dense in the entire complex q-plane. 

Indeed, it suffices to consider the union of the two families 9^ and 6^ + K 2 (the 
latter is the join of 0( s ' p ) with the complete graph on two vertices) and to recall that 
Pg+kM) = q(q-l)---(q-n + l) P G (q - n). 

It is an open question whether the chromatic roots of planar (or perhaps even 
series-parallel) graphs are dense in the disc \q — 1| < 1. I have some partial results 
on this question, but since they are not yet definitive, I shall report them elsewhere. 

The methods of this paper actually prove a result stronger than Theorem 11.21 
namely: 

Theorem 1.4 (a) Fix a complex number vq ^ 0. Then, for eache > and R < oo, 
there exist s G N and a map po: N fl [s , oo) — » N such that for all s > s and 
V > Po( s ), the zeros of Z e ( S , P )( ■ , t>o) come within e of every point in the region 
{q G C: \q + v \ > \v \ and \q\ < R}. 

(b) Fix a complex number q$. Then, for each e > and R < oo, there exist sq G N 
and a map p : N fl [s , oo) — > N such that for all s > s and p > Po(s), the 
zeros of Z e ( s ,p)(qo, ■ ) come within e of every point in the region {v G C: \v\ < 
\qo + v\ and \v\ < R}. 

I thank Roberto Fernandez for posing the question of whether something like Theo- 
rem ED might be true. 

1.3 Sketch of the proof 

The intuition behind Theorem 11.21 is based on recalling the rules for parallel and 
series combination of Potts edges (see Section |21 for details): 

Parallel: v e jj = V\ + t>2 + v \V% (mnemonic: 1 + v multiplies) 

Series: v e g = v \V2/{q + v i + v 2 ) (mnemonic: v /(q + v ) multiplies) 
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In particular, if < \v/ (q + v)\ < 1, then putting a large number s of edges in series 
drives the effective coupling v e g to a small (but nonzero) number; moreover, by small 
perturbations of v and/or q we can give v e jj any phase we please. But then, by putting 
a large number p of such chains in parallel, we can make the resulting v e ff lie anywhere 
in the complex plane we please. In particular, we can make v e jj equal to —q, which 
causes the partition function Zq( s , P ) to be zero. 

To convert this intuition into a proof, I employ a complex-variables result due to 
Beraha, Kahane and Weiss [Ul |H1 E] , which I slightly generalize as follows: Let D 
be a domain (connected open set) in C, and let a%, . . . , a m , . . . , (3 m (m > 2) be 
analytic functions on D, none of which is identically zero. For each integer n > 0, 
define 

m 

f n (z) = J2 a k( z )^(z) n . (1.4) 

k=l 

We are interested in the zero sets 

Z(fn) = {^D: f n (z) = 0} (1.5) 

and in particular in their limit sets as n — > oo: 

liminf Z(f n ) = {z G D: every neighborhood U 3 z has a nonempty intersection 

with all but finitely many of the sets Z(f n )} (1.6) 

limsup i?(/ n ) = {z G D: every neighborhood U 3 z has a nonempty intersection 

with infinitely many of the sets Z(f n )} (1-7) 

Let us call an index k dominant at z if |At(z)| > |A(^)| for alH (1 < I < m); and let 
us write 

-Dfc = {z £ D: k is dominant at z} . (1-8) 
Then the limiting zero sets can be completely characterized as follows: 

Theorem 1.5 Let D be a domain in C, and let . . . ,a m ,/3i, . . . , j3 m (m > 2) be 
analytic functions on D, none of which is identically zero. Let us further assume 
a u no- degenerate- dominance" condition: there do not exist indices k ^ k! such that 
j3 k = ojf3k> for some constant u with \u\ = 1 and such that Dk (= ) has nonempty 
interior. For each integer n > 0, define f n by 

rn 

fc=i 

Then liminf Z(f n ) = limsup Z(f n ), and a point z lies in this set if and only if either 
(a) There is a unique dominant index k at z, and a^{z) = 0; or 
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(b) There are two or more dominant indices at z. 



Note that case (a) consists of isolated points in D, while case (b) consists of curves 
(plus possibly isolated points where all the (3k vanish simultaneously). Beraha- 
Kahane- Weiss considered the special case of Theorem 11.51 in which the /„ are polyno- 
mials satisfying a linear finite-order recurrence relation, and they assumed a slightly 
stronger nondegeneracy condition. Their theorem is all we really need to prove The- 
orems 11.21 and 11.41 but the general result is more natural and its proof is no more 
difficult. Indeed, my proof (see Section |3J) is quite a bit simpler than the original proof 
of Beraha-Kahane- Weiss |7] (though also less powerful in that it gives no information 
on the rate of convergence of the zeros of /„ to their limiting set). 

The next step is to notice that the dichromatic polynomial for the generalized 
theta graph Q^ s ^ has precisely the form (jl.4|) with m = 2: 



Z eis , P )(q,v) 



[{q + v) s + {q- l)v s } p + {q - l)[(q + v) s - v s f 



■jp-i 



;i.9) 



(see Section|21for the easy calculation). It follows from Theorem ll ,5l that when p — > oo 
at fixed s, the zeros of Z e ( a , P ) accumulate where 8 



(9-1) 



q + v 



1 



q + v 



;i.io) 



I then use the following lemma to handle the limit s — > oo: 



Lemma 1.6 Let F x , F 2 ,G be analytic functions on a disc \z\ < R satisfying |G(0)| < 
1 and G ^ constant. Then, for each e > 0, there exists Sq < oo such that for all 
integers s > So the equation 

|l + FiOz)G(z)'| = \l + F 2 (z)G(z) s \ (1.11) 

has a solution in the disc \z\ < e. 

Theorems II . 21 and II .41 are an almost immediate consequence (see Section |S]for details). 



1.4 Plan of this paper 

The plan of this paper is as follows: In Sectional I discuss some identities satisfied 
by the Potts-model partition function Zc{q, {^ e }); in particular, I show what happens 

8 For Theorems ll.2f a') and II. 4f a), the no-degenerate-dominance condition of Theorem 1 1.51 is sat- 
isfied whenever vq ^ 0. For Theorems ll.2f b) and 11.4( b). the no-degenerate-dominance condition is 
satisfied whenever go 7^ 0; but if go = 0, then ZQ( a , P )(qo, ■ ) is identically zero and the assertion is 
trivially true. 
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when a 2-rooted graph (G, x, y) is inserted in place of an edge e* in some other graph 
H . As a special case, I obtain the well-known rules for series and parallel combination 
of Potts edges, which allow one to compute in a straightforward way the Potts-model 
partition function Za[q, {v e }) for any series-parallel graph G. Specializing further, I 
derive the formula (jl.9|) for the dichromatic polynomial Zc{q,v) of the generalized 
theta graphs 0( s > p ). In Section |3] I prove some theorems — which seem to me of 
considerable interest in their own right — concerning the limit sets of zeros for certain 
sequences of analytic functions. As a corollary, I obtain a simple proof of Theorem ll.5l 
In Section 0] I prove a strengthened version of Lemma 11.61 In Section I complete 
the proof of Theorems 11.21 and 11.41 In Section 01 digress to study the real chromatic 
roots of the graphs Q( s ' p \ In Section [7| I discuss some variants of the construction 
employed in this paper. In Section |H] I discuss some open questions. 

The two appendices provide further examples of the power of the identities derived 
in Section |21 In Appendix El I give a simple proof of the Brown-Hickman [25] theorem 
on chromatic roots of large subdivisions. In Appendix El I extend Thomassen's pTSI 
construction concerning the chromatic roots of 2-degenerate graphs. 

2 Some Identities for Potts Models 

In this section I discuss some identities for Potts-model partition functions and use 
them to calculate Zc{q, {v e }) for generalized theta graphs. There are two alternative 
approaches to proving such identities: one is to prove the identity directly for complex 
q, using the Fortuin-Kasteleyn representation (jl.2|) : the other is to prove the identity 
first for positive integer q, using the spin representation (jl.lj) . and then to extend it 
to complex q by arguing that two polynomials (or rational functions) that coincide 
at infinitely many points must be equal. The latter approach is perhaps less elegant, 
but it is often simpler or more intuitive. 

2.1 Restricted Potts- model partition functions for 2-rooted 
graphs 

Let G = {V, E) be a finite graph, and let x, y be distinct vertices of G. We define 
G • xy to be the graph in which x and y are contracted to a single vertex. (N.B.: If 
G contains one or more edges xy, then these edges are not deleted, but become loops 
in G • xy.) There is a canonical one-to-one correspondence between the edges of G 
and the edges of G • xy; for simplicity (though by slight abuse of notation) we denote 
an edge of G and the corresponding edge of G • xy by the same letter. In particular, 
we can apply a given set of edge weights {f e } e gB to both G and G • xy. 
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Let us now define 

z { r y \ q ,{v e }) = £ q k(E,) n«« (2- 1 ) 

E' connects x to y 

Z^KlAve}) = E ^ II* ( 2 - 2 ) 

does not connect x to y 

Note that Z^^ y \q, {v e }) and ZQ^ v \q, {v e }) are polynomials in g and {f e }: the 
lowest-order contribution to Zq* (resp. Z^^^) is of order at least g (resp. at least 
g 2 ); the highest-order contribution to Zq*~* is of order at most q\ v \~ dlst ( x >y) j where 
dist(x, y) is the length of the shortest path from x to y using edges having v e ^ (if no 
such path exists, then Zq*~* v ^ is identically zero); and the highest-order contribution 
to Zq^^ is of order exactly (with coefficient 1, coming from the term E' — 0). 
From p. 2)1 we have trivially 

Z G (q,{v e }) = Z^ y \q, K» + zg^foM) (2-3) 
and almost as trivially 

*Wg, {*.}) = Z^ y \q, {v e }) + q-'Z^iq, {v e }) . (2.4) 

Remark. Let G + xy denote the graph G with an extra edge xy added. Then it is 
not hard to see that (jl.2|) also implies 



J G+xy j ^xy 



- (l+v xy )Z^ y \q,{v e }) + (l + V -^jZ^ y \q,{v e }), (2.5) 

which equals Zq + v xy Zc» X y in agreement with the deletion-contraction formula. In 
particular, when v xy = — 1 we have 

Z G+xy (q, R}, v xy = -1) = ^— - Z ( ^ y \q, {v e }) , (2.6) 

which equals Zq — Zc» X y 

Let now q be an integer > 1, and define the restricted partition function 

Z G,x,y(q, {Ve}', &x, C-y) = £ II [l + ^(^i(e)),^(x 2 (e)))] (2.7) 

o-(x) = o-^ 
cr(y) = o-y 

where a x , a y G {1, . . . , q} and the sum runs over all maps a: V — > {1, . . . , q} satisfying 
a(x) = a x and a(y) = a y . We then have the following refinement of the Fortuin- 
Kasteleyn identity (jl.2j) : 
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Proposition 2.1 

Z G ,x,y(q,{ V e};Vx,(Ty) = A G;Xj y(q, {v e }) + B G>x>y (q, {v e }) 6(<T X , <T y ) (2.8) 

where 

A GiX , y (q,{v e }) = q~ 2 Z^ v \q,{v e }) (2.9a) 

B GiX , y (q,{v e }) = q- l Z^ y \q,{v e }) (2.9b) 

are polynomials in q and {v e }, whose degrees in q are 

degA = \V\ - 2 (2.10a) 
degB < \V\ - 1 - dist(x,y) (2.10b) 

where dist(x, y) is the length of the shortest path from x to y using edges having v e ^ 
(if no such path exists, then B = 0). Moreover, \2.9a /b) define the unique functions 
A G}X ,y and B GjXjV that are polynomials in q and satisfy A2.fy) . 

First Proof. In f)2.7jl . expand out the product over e e E, and let E' C E be 
the set of edges for which the term v e 8(a(xi(e)), a(xz{e))) is taken. Now perform the 
sum over configurations {o~(z)} z ev\{x,y}'- in each connected component of the sub- 
graph (V, E') the spin value o~(z) must be constant. In particular, in each component 
containing x and/or y, the spins must all equal the specified value a x and/or o~ y ; in 
all other components, the spin value is free. Therefore, 



q- 2 Z^ y \q,{v e }) \iOx + o- y 



^G,x,y (q,{v e };a x ,a y 



q- 2 Z ( t y \q, R}) + q-'Z^iQ, if a x = a 



y 



(2.11) 

Finally, the numbers Ac, x ,y and B G x y in (12. 8|) are uniquely defined for each integer 
q > 2; and any polynomials that coincide with (|2.9b /b) at all integers q > 2 must 
coincide everywhere. I 



Second Proof. For each integer q, the S q permutation symmetry of the Potts 
model implies that 

Z G , x , y {q, {v e }; a x , o- y ) = A + BS(a x ,a y ) (2.12) 

for some numbers A and B depending on G, x, y, q and {v e }; moreover, the numbers 
A and B are obviously unique when q > 2. But then, summing ()2.12j) over a x ,a y 
without and with the constraint o x — a y , we get 

Z G = q 2 A + qB (2.13a) 
Z G .xy = qA + qB (2.13b) 
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Hence 



A = Z ° / ZG ' xy = q~ 2 Z^ y) (2.14a) 
q(q-l) 



B = qZ °' xy Z ° = q- l Z { *" y) (2.14b) 
by virtue of (Q and (Q- ■ 

Remark. Extensions of Proposition 12.11 to fc-rooted graphs (for any k > 2) can also 
be derived |14()^ I7H] . I thank Fred Wu for bringing reference |14()j to my attention. 

Let us now consider inserting the 2-rooted graph (G, x, y) in place of an edge e* in 
some other graph H, and let us call the resulting graph H G . We can trivially rewrite 
(M/flU) as 



Zc,x, y {q, {v e }; @xi @y ) = A GiXi y(q, {v e }) 



, B GtX<y (q, {v e }) 

1 + ~a 1 s tt \ x i v) 

A G) x,y{q, \V e \) 



f 2.15a) 



A G ,x, y (q, {v e }) [1 + v efftG>X: y(q, {v e })5(a x , a y )\ (2.15b) 



where 



is a rational function of q and {v e }. We then see that (G,x,y) acts within the graph 
H as a single edge with effective weight t> provided that the partition function Z H 
is multiplied by an overall prefactor Ac lZ , y (q, {v e })- More precisely, 

Z HG{q, {v e }ee(E(H)\e,)UE(G)) = Ao, x ,y (<?, {v e }e£E{G)) Z H (<?, {t>e}eeE(ff)\e» , V e ff) , 

(2.17) 

where v e g = v^c^yiq, {^e}ee-E(G)) replaces v et as an argument of Zh- This follows 
from Proposition ^, ll whenever q is an integer > 1; and the corresponding identity then 
holds for all q, because both sides are rational functions of q that agree at infinitely 
many points. It is also worth noting that the "transmissivity" t e g = v e gj (q + v e g) is 
given by the simple formula 

z<r%,{M) r218 , 

t eff = ~ • \ ZA %) 

The most general version of this construction appears to be the following |137j : 
Let H = (V, E) be a finite undirected graph, and let if be a directed graph obtained 
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by assigning an orientation to each edge of H . For each edge e G E, let G e = 
(V e , E e , x e , y e ) be a 2-rooted finite undirected graph (so that x e , y e G V e with x e ^ y e ) 
equipped with edge weights {fg}g e B e . We denote by G the family {G e } eg £, and we 
denote by H G the undirected graph obtained from H by replacing each edge e G E 
with a copy of the corresponding graph G e , attaching x e to the tail of e and y e to the 
head. Its edge set is thus E = [j E e (disjoint union). We then have: 

Proposition 2.2 Let H = (V,E) and {G e } e( zE be as above. Suppose that 

z G e ,x e , ye (QA v e}eeE e ;cr Xe ,(T ye ) = A G ^ ye {q, {ve}eeE e ) + B Ge!Xe! y e (q, {ve}eeE e ) S(a Xe , a Ve ) , 

(2.19) 

and define 

_ B GejXe! y e (q, {Ve}eeE e ) /„ „m 

AG e ,x e ,y e (q,{Ve}eeEj 

Then 

Z SG (q,{ve}eeE) = I Y[AG e ^ Va (q, {^}e 6 B e ) J x Z H (q, {v eff)e } eeE ) ■ (2.21) 

\eeE J 

In particular, Z^ G does not depend on the orientations of the edges of H . 

A special case of this construction is the subdivision of edges jHSl ESI IH21 EH] , discussed 
further in Appendix 1X1 

2.2 Parallel and series connection 

Important special cases of Proposition 12.11 concern parallel and series connection. 

The case of parallel connection is almost trivial: Let G consist of edges ei, . . . ,e n 
(with corresponding weights V\, . . . , v n ) in parallel between the same pair of vertices 
x, y. Then 



Y[ [1 + Vi5(a x , o-y)\ = 1 + 



i=l 



i=l 



S(a x ,a y ), (2.22) 



so that the parallel edges are equivalent to a single edge with effective weight 

n 

ve ff = n( i+ ^) - 1 ( 2 - 23 ) 

i=l 

or equivalently 

n 
i=l 

13 



That is, for parallel connection, "1 + v multiplies" . 

The case of series connection is a bit less trivial: Let G be a path P n consisting 
of edges ei, . . . , e n (with weights v i, . . . , v n ) between endvertices and n. Then 



£ n^ + ^O 7 *-!' *)] = A + B6(<r ,<r n ) 



(2.25) 



with 



<Tl,...,<T n _l = l 1=1 



A 



!!(?+«*) - n 



i=l 



5 



n 



i=i 



(2.26a) 



(2.26b) 



This formula can easily be proven by induction on n; or it can alternatively be derived 
from (|2.14|) by remembering that 



Z G = Z Pn = qY[(q + Vi) 

i=i 

n n 

Z G .0n = Z Cn = Y[(q + Vi) + (q-l)H 



(2.27) 



(2.28) 



i=i 



i=i 



(formulae that can themselves be proven by induction). In particular, the series edges 
are equivalent to a single edge with effective weight 



V eff = 



B 
A 



qU v i 

i=i 



n n 
i=l i=l 



(2.29) 



or equivalently 



J eff 



q + V e ff 



n 



\q + Vi 



(2.30) 



That is, for series connection, u v/(q + v) multiplies". But one should not forget the 
overall prefactor A given by (|2.26ajl . 

Remarks. 1. Note that applying series reduction (J2.29j) / (|2.3Uj) to a chromatic poly- 
nomial (all Vi = —1) leads to edge weights v e g that are in general unequal and ^ — 1. 
This is further evidence of the value of studying the general Potts-model partition 
function Z G (q, {v e }) with not-necessarily-equal edge weights {v e }. In particular, the 



14 



series-reduction formulae allow an easy proof of the Brown-Hickman theorem on 
chromatic roots of large subdivisions [23]: see Appendix 1X1 

2. The foregoing formulae form the basis for a linear-time [i.e. 0(|V| + \E\)] 
algorithm for computing the Potts- model partition function Z G (q, {v e }) for a series- 
parallel graph G at any fixed set of numbers q and {v e }, and hence a quadratic-time 
[O ( | | 2 + |V||i£|)] algorithm for computing Zc{q, {v e }) as a polynomial in q for fixed 
{v e }. 9 (One needs, as a subroutine, a linear-time algorithm for recognizing whether a 
graph is series-parallel and, if it is, computing a decomposition tree of series-parallel 
reductions: see |127| Sections 2.3 and 3.3].) This is essentially the approach used, 
for the special case of the chromatic polynomial, in [21). More generally, Noble [73] 
has shown how to compute, in linear time, the dichromatic polynomial Zc(q,v) at 
any point (q, v) for graphs of bounded tree-width (the case of series-parallel graphs 
corresponds to tree- width 2); this algorithm presumably extends to handle the general 
Potts-model partition function at any point (q, {v e }). For general graphs, by contrast 
- or even for general planar graphs of maximum degree 4 — even the problem of 
determining whether G is 3-colorable [i.e. whether Pg(3) ^ 0] is NP-complete [J7J 
p. 191]. 



2.3 Potts-model partition function of generalized theta graphs 

Consider a generalized theta graph Sll ... )S consisting of endvertices x, y con- 
nected by p internally disjoint paths of lengths si, . . . , s p > 1. Label the edges 
(1 < i < p, 1 < j < Si) and let Vy be the corresponding weights. Applying the 
series-connection identity (J2.25|) / (|2.26|) on each of the p paths, and then applying the 
parallel-connection identity (|2.22jh we obtain the restricted partition function 

Ze su ..., Sp ,x,y(q,{ V ij} : , (7 x,Vy) = A + B5(a x ,a y ) (2.31) 

with 

(2.32a) 



.4 



n 



ti(9+ - n 



.7=1 



.7 = 1 



B 



n 



.i=i b=i 



lliq + vij) + (g-i)jl 



i=i 



n 



*=i b=i 



Y[(q + Vij) - H 



3=1 



(2.32b) 



In particular, summing ()2.31|) over a x and a y without constraint, we obtain Zq si Bp 
q 2 A + qB: 



9 This is in a computational model where the elementary arithmetic operations (+, — , x, 4-) are 
assumed to take a time of order 1 irrespective of the size of their arguments. 
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Proposition 2.3 The Potts-model partition function for the generalized theta graph 
©si,..., Sp with edge weights {%}i<i< P ,i<j< Sj is 



Ze si ,..., sp (q,{ v n}) 



-(p-1) 



n 



II(g + ««) + (ff-i)il 



i=i 



+ (g-i)II 



*=i u=i 



(2.33) 



In particular, when Vij = v for all i,j, we have the dichromatic polynomial 



Ze sl ,..., 3p (q, v) = \t[[(q + v)« + (q ~ l)v Si ] + (q - 1) [{q 



v \Si _ V H 



. i=l 



i=l 



(2.34) 

And when also s\ = . . . = s p = s, we have 

Z eM {q,v) = q -(p-^{[( q + v y + ( q -iyf + + _„•]*} . (2.35) 

Remarks. 1. Here is an alternate derivation of (j2.31|) / ()2.32|) and ()2.33|) : If in aii „. )Sp 
we contract x to y, we obtain p cycles of lengths Sj joined at a single vertex, so that 

p 

z e Sl ,..., Sp .xy(q, {vij}i<i<p,i<j< Si ) = <T (p_1) Y[Z CH (q, {^}i<j< Si ) 



-(p-i) 



i=i 
p 

n 

i=i 



+ + (g-i)II 



j=i 



(2.36) 



If, on the other hand, we add to Q Sl ,..., Sp an extra edge xy with = —1, we obtain 
p cycles of lengths s» + 1 joined along a single t> = — 1 edge, so that 

^e sl ,..., Sp +a:y(?, {^y}l<i<p,l<j<s<j ^i/ = — 1) 

P 



i=i 



[?(?-i)r M n 



i=i 
p 



( g -l)jJ(g + ^.) - (g-l)JJ 



i=i 



-(p-i) 



(9-1)11 



ii(g+%) - n 



(2.37) 



i=i 



J=i 
16 



j'=i 



The addition-contraction formula Zq = Zc, X y + ZG+xy(v xy = — 1) then gives ()2.33|) ; 
and using (EHl) we obtain (j232>. 

2. In Appendix [B] I use formula (|2.32|) to give an extension (and simple proof) 
of Thomassen's |115t Theorem 3.9] construction concerning the chromatic roots of 
2-degenerate graphs. 

3 Limit Sets of Zeros for Certain Sequences of An- 
alytic Functions 

Let (/„) be a sequence of analytic functions on a domain D C C, and let (a n ) 
be a sequence of positive real numbers such that (|/n| a ™) are uniformly bounded on 
compact subsets of D. We are interested in the zero sets Z(f n ) and in particular in 
their limit sets as n — > oo. We shall relate these sets to the existence and behavior of 
the limit 

u(z) = lim a n log|/ n (z)| . (3.1) 

n — >oo 

We begin with a warm-up result: 

Theorem 3.1 Let D be a domain in C, and let xq G D D M. Let (f n ) be analytic 
functions on D, and let (a n ) be positive real constants such that {\f n \ an ) are uniformly 
bounded on compact subsets of D. Assume further that 

(a) For each x G D D M, f n ( x ) is rea ^ an d > 0. 

(b) For each x G D D M, lim a n log/ n (x) = u{x) exists and is finite. 

n~*oo 

(c) u is not real-analytic at x . 

Then liminf Z(f n ) 3 xq: that is, for all n sufficiently large, there exist zeros z* n of f n 
such that lim z* n = Xq. 

n^oo 

Theorem 13.11 is well-known to workers in mathematical statistical mechanics: it 
is the contrapositive of an observation going back to Yang and Lee |141j concerning 
the genesis of phase transitions (see e.g. |5TH Theorem 4.1, p. 51]). We give its proof 
later in this section. 

The main result of this section, Theorem 13.21 is similar in nature to Theorem 13. II 
the difference is that, by taking the logarithm only of the absolute value of f n (z), we 
need not worry a priori whether f n (z) ^ (since logO = — oo is a legitimate value), 
nor need we worry about potential ambiguities in the branch of the logarithm. 
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Theorem 3.2 Let D be a domain in C, and let zq G D. Let (f n ) be analytic func- 
tions on D, and let (a n ) be positive real constants such that {\f n \ an ) ore uniformly 
bounded on compact subsets of D. Suppose that there does not exist a neighbor- 
hood U 3 z and a function v on U that is either harmonic or else identically — oo 
such that liminf a n log |/ n (^)| < v(z) < limsup a n log |/ n (^)| for all z G U. Then 

n >oo n — >OQ 

z G liminf 2(f n ): that is, for all n sufficiently large, there exist zeros z* of f n such 
that lim z* = z . 

n— >oo 

The proofs of Theorems 13.11 and 13.21 depend crucially on the concept of a normal 
family of analytic functions [721 EH] : 

Definition 3.3 A family JF of analytic functions on a domain D C C is said to 
be normal (in D) if, for every sequence (f n ) C T , there exists a subsequence (f nk ) 
that converges, uniformly on compact subsets of D, either to a (finite-valued) analytic 
function or else to the constant function oo. 

(Some authors omit the possibility of convergence to oo from the definition of normal- 
ity. But the foregoing definition is more useful for our purposes.) An easy covering 
and diagonalization argument |93l Theorem 2.1.2] shows that normality is a local 
property: T is normal in D if and only if, for each z G D, there exists an open disc 
U 3 z in which T is normal. 

The first key result concerning normal families is Montel's (1907) theorem (§51 
pp. 35-36]: if the family T is uniformly bounded on compact subsets of D, then it 
is normal. (In this case, of course, convergence to oo is impossible.) We will need a 
slight extension of this result [95, Example 2.3.9]. For each set A C C, let T^a be 
the family of analytic functions whose values avoid A: 

T iA = {/ analytic on D: f[D] n A = 0} . (3.2) 

We then have: 

Proposition 3.4 If A C C has nonempty interior, then T^a is normal. 

Remarks. 1. Since normality is a local property, this result can be strengthened 
as follows: for a family T to be normal, it suffices that for each z G D there exist an 
open disc U 3 z and a set A v C C with nonempty interior, such that f[U) fl A v = 
for all / G T. 

2. A result vastly stronger than Proposition 13.41 is true: Montel's Critere fonda- 
mental (1912) states that T^a is normal as soon as A contains at least two points. 
Detailed accounts of this deep result can be found in [72^ 19*3]. But we won't need it. 

Proof of Proposition 13.41 Suppose that A contains the disc {w. \w — a\ < e}. 
For each / G T^a-, define f(z) = l/[f(z) — a]. Then the family T = {/: / G J-<£a} 
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is uniformly bounded by 1/e, hence normal by Montel's theorem. So let (f n ) C T^a-, 
and consider the corresponding sequence (/„) C T . By normality of J 7 , there exists a 
subsequence (f nh ) that converges, uniformly on compact subsets of D, to an analytic 
function g. But since the f Uk are nonvanishing on D, Hurwitz's theorem |%T1 p. 262] 
tells us that g is either identically zero or else nonvanishing. It is straightforward 
to show that, in the former case, the corresponding subsequence (f nk ) converges, 
uniformly on compact subsets of D, to the constant function oo; and that in the latter 
case, (fn k ) converges, uniformly on compact subsets of D, to the analytic function 
f — a + 1/g. (For each compact K C D, we have inf \g{z)\ = 5 > 0, so for all 

sufficiently large k we have |/ n J > 5/2 everywhere on K. Now apply the uniform 
continuity of the function w w l/t« on the set where \w\ > 5/2.) | 

We also need the following well-known lemma, which expresses the main idea 
underlying the Vitali-Porter theorem: 

Lemma 3.5 Let T be a normal family of analytic functions on a domain D G C, 
and suppose that the sequence (f n ) C T converges (pointwise either to a complex 
number or to oo) on a set S C D having at least one accumulation point in D. Then 
(f n ) converges, uniformly on compact subsets of D, either to a (finite-valued) analytic 
function or else to the constant function oo. 

Proof. Define f(z) = lim^oo f n (z) for z e S. Now, for any subsequence S = (f nk ), 
there exists (by normality) a subsubsequence S' = (f nk ) that converges, uniformly 
on compact subsets of D, either to a (finite- valued) analytic function gg> or else to 
the constant function g$i = oo. Moreover, we must have gs> \ S = f for all S'. 
Therefore, we either have / = oo, in which case gs> = oo for all S'; or else / is 
everywhere finite-valued, in which case the gs> are all equal to the same analytic 
function g (since S is a determining set for analytic functions on D). It follows that 
(fn) converges, uniformly on compact subsets of D, either to the constant function 
oo (in the first case) or to the analytic function g (in the second). ■ 

We are now ready to prove Theorem 1.111 

Proof of Theorem IH.lt Suppose the contrary, i.e. suppose that there exists an 
e > and an infinite sequence n± < n2 < . . . such that none of the functions f nk has 
a zero in the set D e = {z e C: \z — xq\ < e}; let us also take e small enough so that 
D e C D. Then, since D e is simply connected, log/„ fe is analytic in D e (we take the 
branch that is real on D t fl R), and = a„ fc log/ nfc satisfies lim^oo Uk(x) = u(x) 
for x G D e fl R. Moreover, by hypothesis the functions Re Uk are uniformly bounded 
above on D e , so it follows from Proposition IH.4I that the (uk) are a normal family on 
D e . Lemma 13.51 then implies that the sequence (uk) converges on D e to an analytic 
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function u that extends u. But this contradicts the hypothesis that u is not real- 
analytic at Xq. ■ 

Let us next recall that a real- valued function u on a domain D C C ~ IR 2 is called 
harmonic if it is twice continuously differentiable and satisfies Laplace's equation 

d 2 u d 2 u 

where we write z = x + iy. The key facts (see e.g. [2HJ PP- 122-124]) are: the 
real part of every analytic function is harmonic; and conversely, if the domain D is 
simply connected, then every harmonic function on D is the real part of some analytic 
function. 

One can define normality for families of harmonic functions Section 5.4]: 

Definition 3.6 A family 7i of harmonic functions on a domain D C C is said to be 
normal (in D) if, for every sequence (u n ) C 7i, there exists a subsequence (u nk ) that 
converges, uniformly on compact subsets of D, either to a (finite-valued) harmonic 
function or else to the constant function +oo or — oo. 

A covering and diagonalization argument analogous to Theorem 2.1.2] shows that 
normality is a local property: 7i is normal in D if and only if, for each z G D, there 
exists an open disc U 3 z in which 7i is normal. 

There is a sufficient condition for normality analogous to Proposition 13.41 

Proposition 3.7 Let 7i be a family of harmonic functions on a domain D, which 
are uniformly bounded above on compact subsets of D. Then 7i is normal: that is, for 
every sequence (u n ) C 7i, there exists a subsequence (u nk ) that converges, uniformly 
on compact subsets of D, either to a (finite-valued) harmonic function or else to the 
constant function — oo. 

Proof. Since normality is a local property, it suffices to prove the proposition when 
D is an open disc U. (The only property of U we will really use is its simple con- 
nectedness.) So let (u n ) be a sequence of harmonic functions in U that is uniformly 
bounded above on compact subsets of U . Choose analytic functions /„ on U such 
that Re/ n = u n , and let F n = exp(/ n ). The functions F n are nonvanishing, and 
they are uniformly bounded on compact subsets of U. By Montel's (1907) theorem, 
there exists a subsequence (F nk ) that converges, uniformly on compact subsets of 
U, to an analytic function F. By Hurwitz's theorem, F is either identically zero or 
else nonvanishing. In the former = RelogF nfc = log|F n J tends to — oo 

uniformly on compact subsets of U. In the latter case, u nk tends to the harmonic 
function = Re log F = log \F\ uniformly on compact subsets of U. (For each 
compact K C U, we have inf |.F(/Z)| = 5 > and sup |-F(z)| = M < oo, so for all 
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sufficiently large k we have 5/2 < \F nk \ < 2M everywhere on K. Now apply the 
uniform continuity of the log function on [5/2, 2M].) | 

Remarks. 1. For a slightly different proof, see jH3 Theorems 5.4.2 and 5.4.3]. 

2. Montel's proof [721 Section 23] is not valid: if a subsequence of (f n = u n + zu n ) 
tends to oo, it does not follow that the corresponding subsequence of (\u n \) also tends 
to oo; it could be that (\v n \) does so instead (see e.g. (H3 P- 184]). 

Proof of Theorem 13.21 Suppose the contrary, i.e. suppose that there exists an 
e > and an infinite sequence n\ < n 2 < . . . such that none of the functions f Uk has a 
zero in the set D e = {z G C: \z — z \ < e} C D. Then each function w fc = a nk log \ f nk \ 
is harmonic on D e ; by hypothesis, the functions Uk are uniformly bounded above on 
compact subsets of D t ; so it follows from Proposition 13.71 that the are a normal 
family on D e . Therefore, there exists a subsequence (u^) that converges, uniformly on 
compact subsets of D e , to a function v that is either harmonic on D t or else identically 
— oo. But this contradicts the hypothesis of the theorem. I 



Remark. An argument closely related to that of Theorems 13.11 and 13.21 was used 
previously, in a special context, by Borwein, Chen and Dilcher [201 PP- 79 and 82]. 
I thank Karl Dilcher for informing me of this article. 

We now proceed to the proof of Theorem 11.51 We have already defined 

Dk = {z: k is dominant at z} . (3-4) 
Let us further define the "good sets" 

Vk = {z: k is the unique dominant index at z, and ctk(z) ^ 0} (3.5) 
and the "bad sets" 

Wk = {z: k is the unique dominant index at z, and a^{z) = 0} (3.6) 

X = {z: ^(z) = (3 2 (z) = . . . = f3 m (z) = 0} (3.7) 

Ym = {z: k and I (k ^ I) are dominant indices at z (there may be others) 

and \/3 h (z)\ > 0} (3.8) 

Clearly all the sets (|3.5)) - (j3.8)) are disjoint, except that some of the Y^i may overlap 

m 

at points where there are three or more dominant indices. We shall write V — [J V^, 

fc=i 



21 



m 

W — (J W k and Y — |J Y" fci ; we obviously have the disjoint decomposition 

fe=l l<k<l<m 

D = VUWUXUY. Theorem 11.51 amounts to the statement that 

liminfZ(/„) = limsupZ(/ n ) = W U X U Y . (3.9) 
We begin by collecting some simple facts about these sets: 

Lemma 3.8 Under the hypotheses of Theorem \l.b\ 

(a) Each D k is closed in D. In particular, each dD k has empty interior. 



(b) D k = Vfc U W k U X U I |J Yf-i 1 , and the latter four sets are disjoint. 

(c) For k 7^ I, Dk n Di = X LlY^. In particular, each set X U Yjy is closed. 

(d) For k ^ I, Yki C dDk H dD[. In particular, each Yki has empty interior. 

(e) Each Vk is open. 

(f) Each VkUWk is open. Moreover, V k UW k C D° k C V k UW k UX, where ° denotes 
interior. 

(g) W has no limit points in D. Moreover, for each z G W k there exists e > such 
that the punctured neighborhood {z'\ < \z' — z\ < e} is contained in V k . 

(h) X has no limit points in D. 

(i) W U X UY is a closed set with empty interior, so that V is a dense open subset 
ofD. 

m 

(j) WUXUY= (J 9V k . 

k=l 

(k)Y= (J (YundVtndVi). 

l<k<l<m 

(I) Y is dense-in-itself. 

Proof, (a)-(c) are trivial. 

(d) Let z G Y M C D k n A, so that \(3 k (z ) \ = |A(^o)| > 0. Then fa/ft is analytic 
in a neighborhood U 3 Zq. If fik/Pi equals a constant u, then the no-degenerate- 
dominance condition implies that D k (= D{) must have empty interior, so trivially 
z Q G dD k = dDi. If, on the other hand, p k /Pi is nonconstant, then the open mapping 
theorem implies that Pk/Pi must take values both inside and outside the unit circle 
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in every neighborhood of Zo- The former points belong to and the latter to Df; 
so z G dD k fl dDi. 

(e) is trivial. 

(f) V k U W k is clearly open, hence contained in D k . On the other hand, by (d), 
D° k n F w = for each k, so D° k C 14 U W k U X. 

(g) and (h) are trivial consequences of the assumption that none of the a k or (3 k 
are identically zero. 

(i) We use the fact (valid in arbitrary topological spaces j3H Exercise 1.3.D]) that 
if A is closed, then for every set B we have (A U B)° = (A U B°)°. In particular, if A 
and B have empty interior and at least one of them is closed, then AU B has empty 
interior. Applying this to A = X and B = Y k i, we conclude from (d) and (h) that 
X U Y k \ has empty interior. Moreover, by (c), X U Y ki is closed. Further repeated 
applications, using (g), then lead to the conclusion that W U X U Y is a closed set 
with empty interior. 

(j) By (i), V is open and dense, soWUXUY = V c = dV. On the other hand, 

U Vk ) = U ^^fe f° r an y finite collection of disjoint open sets in any topological 
fc=i / fc=i 
space |HJ Theorem 1 .3.2 (iii) and Exercise 1.3. A]. 

(k) Let z G Y, and let S be the set of indices that are dominant at Zq. Let U be 

an open neighborhood of zo, chosen small enough so that no index in S c is dominant 

at any point of U. By (i), U Pi V is dense in [/. I claim that there must exist at 

least two distinct indices k G S for which U fl V k is nonempty. For suppose that there 

is only one such index, so that U C\V k is dense in U. Then U C D k and z G Y k \ 

for some I ^ k (I & S). Now, shrinking {7 if necessary so that f3 k is nonvanishing in 

U, we have \/3i//3 k \ < 1 in LA, with \(3i(z ) / (3 k (z )\ = 1; so by the maximum modulus 

theorem we must have (3i = w/Jj. for some constant u; with |u;| = 1, contrary to the 

no-degenerate-dominance condition. Hence there are indices k ^ I (k,l G S) such 

that U fl V k and U fl Vi are both nonempty. Since this holds for arbitrarily small 

U 3 z , and since there are finitely many pairs k, I, we must have z G Y k[ fl dV k fl dV\ 

for some pair k ^ I. 

(1) Let zo G F, and let U be any connected open neighborhood of zq, chosen small 
enough so that Un(WUX) = 0; and let U' = U\{z }. If it were true that U'DY = 0, 
then we would have U' C (JfcLi ^fc- But since the \4 are open and disjoint, and U' is 
connected, L/' would have to be contained in one set V k \ but this contradicts the fact 
(k) that every point in Y has nearby points in at least two sets V k and V\. 
■ 

Remarks. 1. It is not necessarily true that Y k i C dV k fl dV\. For example, let 
D — C \ {0} and m = 3, and let /3i(z) = z, ^(z) = 1/z, (3^{z) = 1, ol\{z) = ai{z) = 
a 3 (z) = 1. Then F 12 = F 13 = F 23 = 9Vj = 9V 2 = unit circle, but V 3 = 0. 

2. K cannot have isolated points, but the individual sets Y k \ can: though in a 
neighborhood of zq G Y k i there certainly exist nearby points where \f3 k \ = \/3i\, such 
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points may fail to belong to Y k i because the indices k, I fail to remain dominant. 
Example: fli(z) = 1 + z, ^(z) = 1 — z, /3s(z) = 1. Then Y" 12 = imaginary axis, but 
Y i3 = Y 73 = {0}. 

3. In this paper we won't use parts (j) and (1) of Lemma I3~%1 but they are useful 
facts to know in applications. 

Let us now define 

u n (z) = -\og\f n (z)\ = i Re log /„(*), (3.10) 
n n 

which is well-defined everywhere on D provided that we give it the value — oo at the 
zeros of /„. Clearly, u n is a continuous map from D into MU{— oo}, and is a harmonic 
function on D\Z(f n ). 

We can compute lim u n (z) at nearly every point z G D. Let us define 

n—*oo 

S(z) = {k: a k {z) ^ 0} (3.11) 

P(z) = \ (3.12) 
\ if S(z) =0 

T{z) = {keS(z): \p k (z)\=P(z)} (3.13) 

It then follows easily from the definition of /„ that: 

Lemma 3.9 Under the hypotheses of Theorem M.dc 

(a) For z G V k , lim u n (z) = log |/3fc(^)| > — oo, and the convergence is uniform on 
compact subsets ofV k . 

(b) More generally, if T(z) = {k}, then lim u n (z) = \og\Pk(z)\, though this may 

n—*oo 

equal — oo. 

(c) For z G X, u n (z) = -co for all n, so that lim u n (z) = — oo. More generally, 

n— >oo 

this occurs whenever (3{z) = [i.e. when a k (z)f3 k (z) = for all k]. 

(d) For all z G D, limsupM n (z) < log/?(z). In particular , for z G W k , limsupu n (z) < 

n— >oo n— +oo 

\og\(3 k {z)\. 

Remark. The behavior of the sequence {u n (z)} is subtle at points z where two or 
more dominant indices compete, i.e. where |T(z)| > 2 and f3(z) > 0. But we won't 
need this information. 
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Proof of Theorem 11.51 If z G V k , then choose e > so that D e = {z: \z — z \ < 
e} C V k ; it follows that there exist constants 6 > and M < oo such that 



\otk(z) 


> 


6 


ai(z) 


< 


M 








\Pk(z) 


> 


6 


Pi(z) 
Pk(z) 


< 


1 - 



for all I ^ k 



5 for aU I ^ k 



for all z G D e . Then 



\Uz)\ > \a k (z)\\p k (z)\ n - £|a,(z)||/%(*) 
> <5 n+1 [1 - (m - 1)M(1 - 5) n ] 



(3.14a) 
(3.14b) 
(3.14c) 
(3.14d) 

(3.15a) 
(3.15b) 



for all z G -D e . Therefore, f n is nonvanishing on D e for all sufficiently large n, so that 
z limsupZ(/ n ). 

If zo G we have limsupii n (zo) < log/3(zo) < log lAfcO^o)! by Lemma IQT d). On 

n— >o o 

the other hand, by Lemmas 13.8( g) and 13.9( a) there exists a punctured neighborhood 
of z on which lim u n (z) = log |/3fc(z)|; and this quantity tends to log |/3fc(^o) I as 

n— +oo 

z — ► Zo- Therefore, there does not exist a neighborhood U 3 Zq and a continuous 
function u: [/ — > R U {— oo} satisfying liminf tt n (2) < v(z) < limsup u n (z) for all 

n >oo n — 

z G U. In particular, we can apply Theorem 13 . 21 and conclude that zq G liminf Z(f n ). 

If Zq G X, we have f n (zo) = for all n, so trivially zo £ liminf Z(f n ). 

Finally, if Zq G Y, Lemma l3.8( k) implies that there exist indices k ^ I such that 
zq G Y k i n dVk n dVi. So let C/ be any open neighborhood of z , chosen small enough 
so that (3k and A are nonvanishing on U; then the sets U k = U H V k and Ui = U DVi 
are nonempty. Lemma f3.9f a) implies that 



lim u n (z) 



log \f3 k (z)\ for ^ G U k 
log | AW I iorzeUi 



(3.16) 



U on 



Now there cannot exist a harmonic function v onU that coincides with lim u ny 
U k U Uf. for if there did, then by the uniqueness of harmonic continuation we would 
have IAsWI — IA(^)I — on U, hence (3 k = to>A for some constant uj with |o?| = 1, 
contrary to the no-degenerate-dominance condition. Therefore, by Theorem 13.21 we 
have Zq G liminf Z(f n ). I 
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4 Proof of a Strong Form of Lemma 11.6 



I will prove the following strengthened version of Lemma 11.61 

Theorem 4.1 Let Fi,F2,G be analytic functions on a domain D C C, with G ^ 
constant. For each integer s > 0, define 

U s = {zeD: \l + F 1 (z)G(z) s \ = \l + F 2 (z)G(z) s \} . (4.1) 

Then\immfU s = limsupW, = {z G D: \G{z)\ < 1} \J {z <E D: = \F 2 (z)\} . 

In the proof of Theorem 14.11 I will need the following elementary consequence of 
Taylor's theorem with remainder: 

Lemma 4.2 Let w be a complex number. Then | log(l + w) — w\ < \w\ 2 provided that 
\w\ < 1/4, where "log" denotes the principal branch. 

Proof of Theorem 14. H It is easy to see that limsupW s c {z e D: \G(z)\ < 
1} U {z G D: \F x (z)\ = \F 2 (z)\}: for if z G D with \G(z )\ > 1 and |Fi(zb)| > \F 2 (z )\ 
(say), then we can find a neighborhood U 3 Zq on which \G{z)\ > 1 + 5 and |Fi(z)|- 
\F 2 (z)\ > 5, for some 5 > 0; but then we have |1 + F^G^) 3 ] > |1 + F 2 {z)G{z) s \ for 
all s large enough so that 5(1 + 5) s > 2. 

So it suffices to prove that (a) limmfW, D {z G D: \G(z)\ < 1} and (b) 
liminfW s D{zED: \G{z)\ > 1 and |Fi(z)| = \F 2 (z)\}. 

(a) Assume that zq G D with |G(zo)| < 1- Fix any e > small enough so that 
D e = {z: \z — zq\ < e} C D. Since G is nonconstant, by the open mapping theorem we 
can find Z\ G D e such that |G(^i)| < 1 and G'(zi) ^ 0; and we can then find an open 
set U 3 Z\ with U C U C D € , and a number p < 1, such that G is a homeomorphism 
of [/ onto G[U] G {w. \w\ < p}. If Fi = F 2 , the claim we want to prove is trivial; so 
assume F\ ^ F 2 and choose an open set V C U such that 

(Fi - F 2 )[V] C {w: \w\ > 5 and ipi - ^ < argw < + ^} (4.2) 

for some 5 > and some Since C D e , we have 

M = bu P [|Fi(z)| + |F 2 (z)|] < oo. (4.3) 

Since G[V] is an open subset of {w: \w\ < p}, it contains an arc A r ^ lt o 2 = {w = 
re ie : X < < 6 2 } for some < r < p and Q x < 9 2 . Let C = U fi G^Atf^]; since 
G \ U is a homeomorphism, we know that C is a topological arc (and in particular a 
connected set) contained in V, satisfying G[C] = A r ^ ljd2 . Now take s x = 2n/{6 2 — 9 X ). 
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It follows that for each integer s > s\, the set G S [C] = {G(z) s : z G C} is exactly the 
circle {w: \w\ = r s }. 

Now choose s 2 so that Mp S2 < 1/4. Then, for any integer s > s 2 , the function 

L s (z) = log[l + F 1 (z)G{z) s ] - log[l + F 2 (z)G(z) s ] (4.4) 

(where "log" denotes the principal branch) is analytic on V. Furthermore, by Lemma fOl 
we have 

\L s (z) - [F x {z) - F 2 (z)}G(z) s \ < 2(Mr s ) 2 (4.5) 

for z G C. 

For each integer s > s%, there exist points z£, z~ G C such that G(zf) s = ±r s e~ lipi 
and hence 

±Re{[F 1 (zf)-F 2 (zf)]G(zfY} > ^. (4.6) 
For s > max(si, S2) we then have 

Sv s 

±ReL s (zf) > —= - 2M 2 r 2s . (4.7) 
v2 

Now choose So > max(si,S2) large enough so that 2M 2 r s ° < 5/^/2. It follows that 
for s > so we have ±ReL s (zf) > 0. Since L s is continuous and C is connected, it 
follows that there exists zi°^ G C such that ReL s (zi°^) = 0. This completes the proof 
of assertion (a). 

(b) Assume that z G D with \G(z )\ > 1 and |-Fi(z )| = 1-^2 (^o) I- K ^i(^o) = 
-^2(^0) = 0, the claim is trivial (we have Zq G U s for all s); so assume that |-Fi(z )| = 
I -^2(^0) I > 0. Fix any e > small enough so that D e = {z: \z — z \ < e} C D and 

\G(z)\ > 1 + 5 (4.8a) 
\Fi(z)\ > 5 (4.8b) 
\F 2 (z)\ > 5 (4.8c) 

for all z G D e , for some 5 > 0. We then consider two cases: 

Case 1: Fx = ujF 2 for some uj G C (obviously = 1). We need to show that for 
all s sufficiently large, there exists z G D e such that 

ll + F^z^Giz)' 8 ] = \l + F 2 (z)- l G(z)- s \. (4.9) 

But this follows from part (a) applied to the functions 1/ F 2 and 1/G. 

Case 2: There does not exist uj G C such that Fi = ujF 2 . In this case the function 

f(z) = F 1 {z )F 2 (z) - F 2 {z Q )F x {z) (4.10) 
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has a zero at z = Zq but is not identically vanishing. So it follows from Rouche's 
theorem that for all sufficiently large s, the equation 

F 1 (z )F 2 (z) - F 2 (z )F 1 (z) + [F^zo) - F 2 (z )]G(z)- s = (4.11) 

has a solution z G D e . But this means that 

F 1 (z)+G(z)- S F 2 (z)+G(z)- S 



Ft(zo) F 2 (z ) 

hence that z I 



(4.12) 



5 Completion of the Proof of Theorem 11.41 

We can now put together all the ingredients to complete the proof of Theorem II. 41 
(and hence also of the weaker Theorem I1.2J1 along the lines sketched in Section 11.31 
We shall give the details for part (a); the proof of (b) is almost identical, mutatis 
mutandis. 

Fix a complex number vq 7^ 0, and fix e > and R < 00. Since the open region 
D = {q G C: |(? + i>o| > |^o| and |g| < R} has compact closure D = {q G C: |<7 + i>q| > 
\vq\ and \q\ < R}, we can choose finitely many points q±, . . . , q n G D such that each 
point q G D lies within a distance e/3 of at least one 

Now, for each i (1 < i < n), apply Lemma 1 1.61 (or the stronger Theorem I4.1jl 
with z = q — qi, F\(z) = q — 1, F 2 (z) = —1, G(z) = Vo/(q + v ). [The hypothesis 
v 7^ guarantees that G ^ constant, and the hypothesis qi G D guarantees that 
|G(0)| < 1.] We conclude that, for each i, there exists an integer s 0ti < 00 such that 
for all integers s > s 0ti the equation 

i+(«-i)f-^y = 

[cf. f|1.10|l ] has a solution in the disc \q — < e/3. Let s Q = max s 0i . 

l<i<n ' 

Now fix any s > sq, and let q[, . . . , q' n be solutions of (|5.1jl satisfying \q[ — qi\ < e/3 
and q\ 7^ for all i. Apply Theorem 11.51 (or the original Beraha-Kahane-Weiss 
theorem) with m = 2, n = p, z = q, cei(z) = 1, a 2 (z) = q — 1, (3\(z) — (q + v ) s + 
(q — l)t>0) Pi(z) = (q + vo) s — Vq. [The nondegeneracy condition holds because Vq 7^ 0.] 
We conclude that the zeros of 

<f -%< S , P > (q, v) = [(q + Vo y + (q- l)v&> + (q - 1) [(q + v ) s - (5.2) 

[cf. 1)1.9)1 ] accumulate (in lim inf as well as lim sup sense) at each point of the solution 
set of 1)5.1)1 . In particular, for each i, there exists a zero of ()5.2jl in the disc \q — q[\ < 



'-I — ) 
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min(e/3, |g-|). [Since this disc does not contain 0, the given q must be a zero of 
Z e ( s , p) (q,Vo).] Set p = max p 0i . 

l<i<n 

We have shown that for each p > p , there exist zeros q" , . . . , q" of Zq( S , P ) ( • , i>o) 
such that each point q G D lies within distance e of at least one q". This completes 
the proof. 

6 Real Chromatic Roots of the Graphs 

Though this paper is primarily concerned with the complex chromatic roots of 
the graphs 0( s,p ) [and more generally, with the complex roots of their dichromatic 
polynomials Z e ( a , P ) (q, v)], let me digress to discuss their real chromatic roots. The 
chromatic polynomial P e { S , P )(q) is obtained by specializing ()1.9|) to v — — 1: 



i(q - i) s + (q- i)(-i)T_+ (q - - i) s - (-i)T 

q' j 



■ ^p)(q) = — x • (6-1) 



It is useful to write y = 1 — q, so that 

F(y) = (-lr+Y-^Pe^g) = y(y s - If - (y s - V) v ■ (6.2) 

Recall that 0( s,p ) has n = (s — l)p + 2 vertices. I write P^\q) to denote the j th 
derivative of Pc(q)- The real chromatic roots of the graphs 0( s,p ) are then fully 
characterized by following result: 

Proposition 6.1 Fix integers s,p > 2 and set n = (s — l)p + 2. Then: 

(a) (-l) n+j P^ SiP) (q) > for q < and < j < n. 

(b) P 0Cs , P) (O) = and (-l) n P^ p) (0) = (s- If - < 0. 
f C ; (-lfP e(s „)(g) <0/or0<g<l. 

^ Pe( SlP )(l) = and (-1)^^(1) = 1 > 0. 

(^ej zs even (i.e. s is odd or p is even or both), then P @ { S , P )(q) > for q > 1. 

(e ' ) Ifn is odd (i.e. s is even and p is odd), then there exists qi s ' p ^ £ [1.4301597 . . . , 2) 
such that: 



(e[) P & (s, p) (q) < for 1< q < q 

Dl 



( s >p) . 



(e' 2 ) P 0( s, P) (qi s ' p) ) = and P'(qi s ' p) ) > 0; and 



(e'J P e (s, P ){q) > and P' (s , p) {q) > for q > q 



( s .p) 
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Here 1.4301597 ■ ■ ■ (— qi 2 ' 3 ^ ) is shorthand for the unique real root of q 3 — 5q 2 + 
lOg — 7 = 0. Moreover, 

(e' a ) qi s ' p ^ is a strictly increasing function of s and p; 
(e' h ) lim qi = 2 for each p; and 

(e' c ) lim qi s ' p ^ = qi s '°°\ where gl s ' 00 ' ) G [1,2) is the unique nonzero real root of 

p—>oo 

2(g-l) s + g-2 = 0. 

The signs of Pdq) and its derivatives in (a)-(d) are of course general facts about the 
chromatic polynomial of any 2-connected graph |1331 Theorem 2], but it is amusing 
to verify them explicitly for @( s ' p \ 

Proof, (a) We have 

F(y) = (y - l)(y s - If + [(y s - If - (y s - yf] (6.3a) 

p-i 

( y s -iy + - 1) V - v) 1 ^ 1 " 

k=0 



(y-i) 



(6.3b) 



so clearly F and all its derivatives are nonnegative for y > 1. If we now write 
Pewb) = ELi^V, we see that (-l) n - k a[ s ' p) = F^ 1 ) (1) / (k + p - 1)! > 0; 
and of course g4 = 1 > 0. The claims follow. 

(b) An easy calculation. 

(c) For < y < 1, we have (1 — y s )/{y — y s ) > 1/y > 1 and hence [(1 — y s )/(y — 
y s )f > l/y, so that (-l)PF(y) > 0. 

(d) An easy calculation. 

(e) For y < and p even, we manifestly have F(y) < 0. For y = —z < and s 
odd, we have (— l) p+1 F(y) = z(z s + l)' p + (z s — z) p . For z > 1 this is manifestly > 0; 
for < z < 1 it is > z — z p > as well. 

(e') For z = q — 1 > with s even and p odd, we can write P e ( 3 , P )(g) = (1 + 
z)G{z) p [l -H(z)p) where 

G(z) = (6.4a) 

Now i?(^) is a strictly decreasing function of z on (0,oo), which runs from +oo at 
z = to at z = 1 to — oo at z = +oo; therefore, the equation H(z) = 1 has 
a unique solution zi s ' p ^ 6 (0, oo), which lies in fact in (0,1). Moreover, #(2) is a 
strictly increasing function of s and p when z G (0, 1), so z[ s ^ is a strictly increasing 
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function of s and p; and it is easy to see that it has properties (e'b) and (e' c ). Finally, 
G(z) is a strictly increasing function of z on (0, oo), so -Pq( 3 , p) (<?) > for q > g* . I 

Remarks. 1. For any n- vertex connected graph G and any vertex v of G, the quantity 
Tc(l, 0) = (—l) n+1 Pq(0) counts the acyclic orientations of G in which v is the unique 

source jlHl Theorem 7.3] (23 Proposition 6.3.18] jlH]. Taking v to be one of the 

p p 

endvertices of 6 sl) ... jS „, it is easy to show that this quantity equals FJ — n( s « ~~ -0> 

exactly as computed from the chromatic polynomial. 

2. For any n-vertex connected graph G ^ K 2 and any edge e = (w) of G, the 
quantity [<9T G (x, y)/cfc](0, 0) = (—l) n P' G {l) counts the acyclic orientations of G in 
which u is the unique source and v is the unique sink ^DJ Theorem 7.2] [23 Exercise 
6.35] jlH]. It also equals half the number of totally cyclic orientations of G (i.e. 
orientations in which every edge of G belongs to some directed cycle) in which every 
directed cycle uses e [HU Theorem 8.2] [23 Proposition 6.2.12 and Example 6.3.29]. 
Taking u to be one of the endvertices of Sli ... jSp and v one of its neighbors, it is easy 
to see that there is exactly one orientation of the first kind (and it is acyclic), and 
two of the second. 

3. In part (e) we cannot assert that P e ( S , P )(q) is increasing for q > 1. Indeed, 
for many pairs (s,p) — including (2,4), (2,6), (4,6), (5,6), (5,7), ... — there is an 
interval l<qi<q<q 2 <2 where -Pq (s , p ) ((?) < 0. 

4. In part (e'), P e (s, P )(q) is not necessarily convex on (1,2): for example, for 
(s,p) = (2,7) there are inflection points at q ~ 1.282916,1.405642. However, nu- 
merical calculations for small s,p suggest that P@( s , P )(g) is convex on the interval 
[qi s ' p \ oo). Indeed, it appears that all the derivatives of P e ( S , P )(q) are strictly positive 
at q — qi s ' p \ except that for s = 2 and p > 7 the (n — l) st derivative is < 0. 

Though I am unable to say much about the monotonicity or convexity of Pq( s , p ) (q) 
for 1 < q < 2 (cf. Remarks 3 and 4 above), the situation simplifies dramatically for 
q > 2 for all series-parallel graphs. We begin with the following: 

Definition 6.2 Let P be a degree-n polynomial with real coefficients, and let a e K.. 
We say that P is strongly nonnegative (resp. strongly positive^ at a in case P^\a) > 
(resp. P( j \a) > 0) for all < j < n. 

Note that if P is strongly nonnegative at a and is not identically zero, then P is 
strongly positive on (a, oo); and conversely, if P is strongly nonnegative on (a, oo), 
then it is also strongly nonnegative at a. Note also that if P and Q are strongly 
nonnegative (resp. strongly positive) at a, then so are P + Q and PQ. I suspect 
that strong positivity is the appropriate concept for many (though perhaps not all) 
theorems asserting upper zero- free intervals for chromatic polynomials [T7| I136tlll5j . 
Here, at any rate, is one example: 
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Proposition 6.3 LetG = (V, E) be a finite graph equipped with edge weights {v e } e£ E, m 
let x, y be distinct vertices of G; and let 

) = A GiXiy (q,{v e }) + B G)X)y (q,{v e })6(a x ,cr y ) (6.5) 

be its restricted partition function. Suppose that the two-terminal graph (G, x, y) is 
series-parallel 10 , and that v e > — 1 for all edges e G E. Then A and A + B are 
strongly nonnegative at q = 2. 

Corollary 6.4 Let G be a series-parallel graph equipped with edge weights {v e } e( zE 
satisfying v e > — 1 for all edges e. Then Z G (q, {v e })/q [and hence also Z G (q, {v e })] 
is strongly nonnegative at q = 2. In particular, P G (q) > for q > 2. 

The proof of Proposition 16.31 is an easy application of the formulae ()2.23j) and 
(12.26)1 for parallel and series connection, which can be rewritten as 

({A 1 ,B X ),{A 2 ,B 2 )) pa ^ lcl (A^,B^) = (A 1 A 2 ,A 1 B 2 + A 2 B 1 + B 1 B 2 ) (6.6) 

B,), (A 2 , 5 2 )) se A es (A se \ 5 scr ) ee {A\B 2 + A 2 B X + qA 1 A 2 , B X B 2 ) (6.7) 

Lemma 6.5 If A\, A\ + B\, A 2 and A 2 + B 2 are all strongly nonnegative at q = 2, 
then so are A par , A par + _B par ; A sei and A sei + B sev . 



PROOF. A trivial computation: 



^par _|_ j^par 
A ser 
yjser _|_ J^^er 



A X A 2 (6.8a) 

(A 1 + B 1 )(A 2 + B 2 ) (6.8b) 

Ax(A 2 + B 2 ) + A 2 (Ai + B x ) + (q- 2)A X A 2 (6.8c) 

(A 1 + B 1 )(A 2 + B 2 ) + (q-l)A 1 A 2 (6.8d) 



Proposition 16.31 then follows by induction, starting from the fact that it holds when 
G is a single edge (xy) (for which A = 1 and A + B = l-\-v e >0). Corollary 16.41 for 
2-connected series-parallel graphs G then follows from Proposition 16.31 by choosing 
some pair of vertices between which the graph is series-parallel, and using Z G /q = 

10 Note that this is stronger than requiring simply that G be series-parallel; we require that G 
with the two fixed terminals x and y be transformable to a single edge by a sequence of series and 
parallel reductions. 
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qA + B = (A + B) + (g — 1)A; it then holds for general series-parallel graphs by 
factorization into 2-connected components. 

Remark. Thomassen |115| Theorem 3.4] shows that for any integer k > 2, the set 
Mj. of real chromatic roots of graphs of tree-width < k consists of 0, 1 and a dense 
subset of the interval (32/27, k}. Indeed, Thomassen's argument (which is based on 
|115| Theorem 3.3]) actually proves the stronger result that -Pg(<?)/? [and hence also 
Pg(q)] is strongly nonnegative at q — k. The case k = 2 corresponds to series-parallel 
graphs. 

7 Variants of the Construction 

7.1 Adding an edge to G)( s >p) 

I recently proved that if G is a graph of maximum degree A, then all the chromatic 
roots of G lie in the disc \q\ < 7.963907A jl 1 1| Corollary 5.3 and Proposition 5. 4]. 11 
Moreover, if all but one of the vertices of G have degree < A, then the chromatic 
roots of G are still bounded, namely they lie in the disc \q\ < 7.963907A + 1 [111, 
Corollary 6.4]. The graphs 0( s ' p ) show that no analogous theorem can hold if "all 
but one" is replaced by "all but two" (not even when A = 2 and G is series-parallel), 
for in that case the chromatic roots can spread over the whole complex plane (except 
perhaps the disc \q — 1| < 1). 

There is at least one instance in which the presence of several vertices of large 
degree cannot cause large chromatic roots, namely when those vertices form an N- 
clique within which every pair of vertices is connected by a v = — 1 edge. In that 
case, if the remaining vertices of G have degree < A, then all the chromatic roots of 
G lie in the disc \q\ < 7.963907A + N [TJH Theorem 6.3(b)]. 

The condition here that the clique edges have v = — 1 is crucial. Indeed, consider 
the graph 0( s ' p ) + xy (where x and y are the endvertices of 0( s,p )) in which the edge 
xy is assigned a weight v xy and all the other edges have weight v. The corresponding 
Potts-model partition function is easily calculated from ()2.5|) . ()2.8|) and (j2.31|) / ()2.32|) : 

7 , v (1 + v xy )[(q + v) s + (q- 1KF + (g ~ l)[(g + vY ~ v S ] P 

q (7-1) 

n This result is, in fact, the specialization to chromatic polynomials of a more general bound on the 
zeros of the Potts-model partition function Zo{q, {ve}) throughout the "complex antiferromagnetic 
regime", i.e. when |1 + v e \ < 1 for all edges e. The other two results of (11 1| cited later in this 
subsection are likewise specializations of bounds valid throughout the complex antiferromagnetic 
regime. See for details. 
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which reduces to (I2.35|) when v xy = and to ()2.37|) when v xy = — 1. Provided that 
v xy 7^ — 1, Theorem 11.51 still applies and the zeros accumulate as s,p — > oo in the 
entire region where \v/(q + v)\ < 1. 

7.2 The dual construction 

Let Cn be the cycle C n with each edge replaced by r edges in parallel. Note that 
the dual graph of 0( s > p ) is Cp\ and vice versa. The dichromatic polynomial of Cn^ 
is easily computed from the dichromatic polynomial ()2.28j) of C n together with the 
parallel reduction formula ([2.22)1 : 

Z c (Mq,v) = [ q + (l+v) r -l] n + (q-l)[(l + v y-l] n . (7.2) 

It follows from Theorem ll . 5l that when n — > oo at fixed r, the zeros of Z ( r ) accumulate 
where 

|1 + (g-l)(l+^)- r | = \l-(l+v)- r \. (7.3) 

Here the G ^ constant condition of Lemma fl. 61 means that we cannot fix v and use q 
as the variable. However, we can fix q = q ^ and use v as the variable 12 ; we then 
have: 

Theorem 7.1 Fix complex numbers qo,Vo satisfying q^ ^ and |1 + Vq\ > 1. Then, 
for each e > 0, there exist numbers r < oo and no(r) < oo such that for allr > r and 
n > n (r), the dichromatic polynomial Z ( r )(qo, ■ ) has a zero in the disc \v — v \ < e. 

We can also slice the (q, v )-space in other ways; for example, we can fix w ^ 
and take v = q/u : 

Theorem 7.2 Fix nonzero complex numbers qo,Uo satisfying \uq\ < \qo + Uq\. Then, 
for each e > ; there exist numbers r < oo and no(r) < oo such that for all r > r 
and n > n (r), the polynomial Z ( r )(q, q/u ) has a zero in the disc \q — q \ < e. 

In particular, if\qo—l\ > 1, then the flow polynomial F ( r )(q) = (—l) nr q~ n Z( r )(q,—q) 
has a zero in the disc \q — qo\ < e. 

Since the G ^ constant condition forbids us to use q here as the variable with 
v fixed, these results are somewhat less interesting (at least for direct application 
to chromatic polynomials) than Theorem 11.21 However, they become useful when 
we insert nontrivial 2- rooted graphs in place of the edges of Cn '\ But we defer this 
construction to a separate paper. 

12 Note that the no-degenerate-dominance condition of Theorem ll. 51 is satisfied provided that the 
fixed value of q is nonzero. 
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8 Discussion 



It now appears that there are three quite distinct theories of chromatic roots: 

(a) the theory of integer chromatic roots [i.e. the combinatorial theory of the chro- 
matic number x(G)]; 

(b) the theory of real chromatic roots (i.e. the theory of zero-free and zero-dense 
intervals for various classes of graphs [T71 E£2 EH HSS1 E3 EM EES IH ITTflj ): 

and 

(c) the theory of complex chromatic roots. 13 

Woodall |1331 Theorem 8] has shown that the complete bipartite graphs K ni n2 , in 
the limit n 2 — > oo with n\ fixed, have real chromatic roots arbitrarily close to all the 
integers from 2 through |_fta/2j, even though their only integer chromatic roots are 
of course and 1. More generally, the joins G = K niiU2 + K n , which have chromatic 
number x{G) = n + 2, have real chromatic roots arbitrarily close to all the integers 
from n + 2 through n + |_^i/2_|, even though their only integer chromatic roots are 
0, 1, . . . , n + 1. This suggests that the real roots of the chromatic polynomial have 
nothing much to do with the integer roots. In particular, the real roots of Pg{q) 
cannot be bounded in terms of the chromatic number alone. 

Likewise, the results of the present paper show that the generalized theta graphs 
q(s,p) have complex chromatic roots arbitrarily close to any point in C\{|g— 1| < 1}, 
even though their only real chromatic roots are 0, 1 and q[ s ' p) G [1.4301597 2) 
and their only integer roots are and 1. This suggests that the complex roots of the 
chromatic polynomial have nothing much to do with the real or integer roots. 

In particular, planar graphs (such as B^ s,p - ) ) can have complex chromatic roots 
densely surrounding the point 4. 14 This suggests that it is unlikely that there will 
ever be a "complex-analysis" proof of the Four-Color Theorem: quite simply, q = 4 

13 There is also a (small) fourth theory, namely the "number-theoretic" (or "algebraic") theory 
of chromatic roots, which exploits the fact that the chromatic polynomial is a monic polynomial 

(k) 

with integer coefficients. I include here the proof that the generalized Beraha numbers B„ = 
4cos 2 (fc7r/n) for n = 5,7,8,9 and n > 11, with k coprime to n, are never chromatic roots |92) 
(see |123) for earlier related results); Tutte's golden-ratio theorems for plane triangulations |120l 
121, 123; and speculations concerning the accumulation of chromatic roots at Beraha numbers 
B n = 4cos 2 (7r/n) for certain sequences of planar graphs QUI E2 EH El EH EI EM 123 EH ESI El 
EUE3E3EEIEHlEHlEZIEElE2EllS3j- This theory is, however, as yet rather undeveloped. 

14 The fact that planar graphs can have chromatic roots arbitrarily close to 4 was established 
already two decades ago in the seminal paper of Beraha and Kahane They considered 4 x n 
strips of the triangular lattice with periodic boundary conditions in the "short" direction (and extra 
sites at both ends of the "long" direction, but this is inessential) and found that the chromatic roots 
accumulate as n — > oo on a curve passing through q = 4. 
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behaves very differently from complex values of q arbitrarily close to 4. On the 
other hand, it is conceivable that there may some day be a "real-analysis" proof of 
the Four-Color Theorem in the same sense that Birkhoff and Lewis [T71 pp. 413- 
415] (see also |136| 1115] ) provided a "real-analysis" proof of the Five-Color Theorem: 
namely, it may be possible to prove that Pg(q) > for real q > 4, for all loopless 
planar graphs G. Indeed, Birkhoff and Lewis [T7[ p. 413] have conjectured that 
Pg{q) /[q(q ~ !)(<? ~~ 2)] — (q — 3) n ~ 3 is strongly nonnegative (see Definition 16. 2|) at 
q = 4 for all loopless planar graphs G with n vertices. 

We conclude with a few disconnected remarks: 

1) Let G% denote the class of graphs G such that every sub contraction (= minor) 
of G has a vertex of degree < k. Woodall |136l Question 1] asked whether all the 
(complex) chromatic roots of G e are bounded in modulus by k. For k > 2, the 
answer is no, since 

generalized theta graphs C series-parallel graphs = (8.1) 

and we have shown that the chromatic roots of generalized theta graphs 0( s ' p ) are 
unbounded. (In fact, the family 0( 2 ' p ) ~ K 2iP already has unbounded chromatic roots: 
see HH.) 

2) Read and Royle [HS1 PP- 1027-1028] have observed empirically that for cubic 
graphs with a fixed number of vertices, the graphs with high girth tend to contribute 
the chromatic roots with smallest real part. It would be valuable to obtain a better 
understanding of this phenomenon. A first conjecture might be that there is a uni- 
versal lower bound on the real part of a chromatic root in terms of the girth, but 
this is false: the graphs 0( 3 ' p ) all have girth = circumference = 6, but their chromatic 
roots accumulate as p — > oo on the hyperbola 

(Img) 2 - (Reg - |) 2 = \ (8.2) 

[set s = 3 and v = — 1 in (jl.lOj) ]. which contains points with arbitrarily negative real 
part. 

3) The construction employed in this paper is based on concatenating 2-rooted 
subgraphs (here single edges, but that can be generalized) in order to create a larger 
2-rooted subgraph with some desired value of the effective coupling v e g. The situation 
may thus be very different for 3-connected graphs, in which 2-rooted subgraphs of 
more than a single edge cannot occur. It is thus reasonable to ask: What is the closure 
of the set of chromatic roots of 3-connected graphs? The answer, as Bill Jackson has 
pointed out to me, is "still the whole complex plane". To see this, consider the 
graphs 0( s,p ) + K n for fixed n and varying s,p: they are {n + 2)-connected, and 
their chromatic roots (taken together) are dense in the region \q — (n+ 1)| > 1. In 
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particular, considering both n and n + 2, the chromatic roots are dense in the whole 
complex plane. So even arbitrarily high connectedness does not, by itself, stop the 
chromatic roots from being dense in the whole complex plane. 

On the other hand, the graphs @( s ' p ) + K n are non-planar whenever n > 1 and 
p > 3. This suggests posing a more restricted question: 

Question 8.1 What is the closure of the set of chromatic roots for 3-connected planar 
graphs ? 

Here the answer may well be much smaller than CorC\{|g — 1| < 1}. But it 
will certainly not be a bounded set, for as Bill Jackson has pointed out to me, the 
bipyramids C n + K 2 are 4-connected plane triangulations, but their chromatic roots 
are unbounded jH QUI EU • 

4) The chromatic roots of generalized theta graphs have been further studied in 
|1U5| |2H HE!- In particular, Shrock and Tsai |lU5t Section 3] plot the limiting set 
(jl.lUJl of chromatic roots as p — > oo for fixed s = 3,4, 5. Brown, Hickman, Sokal and 
Wagner prove that the chromatic roots of S1) ... )S are bounded in magnitude by 
[1 + o(l)] pj logp, uniformly in the path lengths S\, . . . , s p ; moreover, they prove that 
for 0( 2 ' p ) ~ i^2 iP this bound is sharp. 
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A Simple Proof of the Brown— Hickman Theorem 
on Chromatic Roots of Large Subdivisions 

In this appendix I use the formula (|2.25|) / ()2.26j) for series reduction of Potts edges 
to give a simple proof of the Brown-Hickman [25] theorem on chromatic roots of large 
subdivisions. This proof gives a further illustration of the utility of considering the 
general Potts-model partition function Z@(q, {v e }) with not-necessarily-equal edge 
weights v e , even if one is ultimately interested in the dichromatic polynomial Z G (q, v) 
or the chromatic polynomial P G {q)- 

Lemma A.l Let G be a graph of n vertices and m edges, and let all the edge weights 
satisfy \v e \ < 5. Then all the roots ofZ G (q, {v e }) lie in the disc \q\ < max(C, C 1 ^ 71 ^) 
where C = (1 + 5) m - 1. 

PROOF. From (jl.2|) we have Z G (q, {v e }) = Y^k=i a kq k with a n = 1 and Ylk=i l afe l — 
(1 + 5) m - 1 = C. Then ^ ' 

^ t { C\q\ n ~ l if \q\ > 1 ) 

fc= i L^kl if \q\ < i J 

So Z G ^ if |g| > max(C, C 1 /^). I 

Remark. This bound is far from being quantitatively optimal; much better bounds, 
which depend on the maximum degree A but not on n or m, can be obtained using 
the methods of But Lemma lA.ll is good enough for our present purposes, as 

we do not seek any uniformity in G. 

Given a finite graph G = (V,E) and a family of integers k = {k e } e( zE > 1, we 
define G k to be the graph in which the edge e of G is subdivided into k e edges in 
series. 

Theorem A.2 Let G be a finite graph, let v be a complex number, and let e > 0. 
Then there exists K < oo such that k e > K for all edges e implies that all the roots 
of Z G k(q, v) lie in the disc \q + v\<{l + e)\v\. 

Proof. Let G have n vertices and m edges. The claim is trivial if v = 0, so assume 
that v 0. By the series reduction formula (|2.25jl / (j2.26j) . we have 

Z(f(q,v) = [II iq + v)ke ~ vke ) x Z G (q,{v e , eff }) (A.2) 

\e£E q J 
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where ^ 

Ve > eff = (q + - ^ • (A ' 3) 



The zeros of the prefactor in (jA.2|) all lie on the circle \q + v \ = |t>|. So let us 
assume that \(q + v)/v\ = R > 1 + e; we shall show that Z G (q, {f e>ej (f}) 7^ provided 
that K is chosen sufficiently large (depending only on e, n, m and |v|). We have 
e < R — 1 < |g/t>| < i? + 1 and hence 

for all edges e [since (R + 1)/(R K — 1) is a decreasing function of R for R > 1]. It 
follows from Lemma I A. II that Zc{q, {v e ,eff}) 7^ provided that we choose K large 
enough so that 5 < (1 + a) 1 / m — 1 where a = min[e|t> |, (e|t> l) 1 ^" -1 ^]. ■ 

The Brown-Hickman theorem is the special case v — — 1. 



B Chromatic Roots of 2-Degenerate Graphs 

Let us recall that the degeneracy number of a graph G is defined by 

D(G) = max5(#) , (B.l) 

where the max runs over all subgraphs H of G, and 6(H) is the minimum degree of 
H. A graph G is said to be k-degenerate if D(G) < k. Equivalently, G is fc-degenerate 
if there exists an ordering x±, X2, ■ ■ ■ , x n of the vertex set of G such that 

'^.••{.n..,-, , ;'•'•.;) < k ( B - 2 ) 

for 1 < j < n, where G \ {x±, X2, . . . ,Xj} denotes the induced subgraph. 

It is an elementary fact jT^l Theorem V.l] that all of the integer chromatic roots 
of G lie in the interval [0, D(G)}, i.e. the chromatic number x{G) is & t most D(G) + 1. 
On the other hand, Thomassen |115t Theorem 3.9] has shown that the real chromatic 
roots cannot be bounded in terms of the degeneracy number: indeed, there are 2- 
degenerate graphs with arbitrarily large real chromatic roots. Here we shall use 
Proposition 12.21 and equation ()2.32|) to give a simple proof of the formula underlying 
(a slight generalization of) Thomassen's construction, and to exhibit its exact range 
of applicability. 

One sees immediately from (J2.16j) / (J2.32j) that the 2-rooted generalized theta graph 
@0>p) (rooted at the endvertices) has an effective coupling 
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where 

Mq) . fa- 1 ';-'- 1 )' (b.4, 

CM) = fa- 1 '^';- 1 "- 1 '' (B.5) 

So let Qs C C be the set on which |1 + i> ej gr(?)| > 1, i.e. 

Q s = {g G C: |(g - l) s + (g - l)(-l) s | > |(g - l) s - (-l) s |} , (B.6) 
so that in particular 

Q 2 = {g: Reg > §} (B.7a) 

Q 3 = {g: (Img) 2 - (Reg-|) 2 > \} (B.7b) 

We then have the following easy extension of Thomassen's result: 

Theorem B.l Let G be a loopless graph, and let s be an integer > 2. Then there 
exist 2-degenerate (hence 3-colorable) graphs {G'(p)}^ =1 and an integer k > such 
that 

(a) lim C s (q)~ kp Pqi t v \{q) = Pg(q) uniformly on compact subsets of Q s . 

p^oo 

(b) If G is connected, then so are the G'{p). 

(c) If s is even and G is bipartite, then so are the G'{p). 

(d) If G is embeddable on a surface S, then so are the G'{p). 
In particular: 

(i) V % E Q s is a zero of Pq, then for all e > and all sufficiently large p 
(depending on e) Pg'(j>) has a zero in the disc \q — q \ < e. 

(ii) If g G Q s fl R is a zero of Pq of odd multiplicity, then for all e > and all 
sufficiently large p (depending one) Pg'( p ) has a zero in the interval (g — e, go+e)- 

Proof. We can assume without loss of generality that G = (V, E) has no multiple 
edges (this is just to simplify the notation). For each x G V, define 

f \x\ if x is an isolated vertex of G 

^ {{x, y): xy G E) otherwise 
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and let V = U x&v V' x - Let E' = {(x,y)(y,x): xy G E}, and let E" be a set of edges 
(with vertex set V) chosen so that the vertex sets of the connected components of the 
graph (V, E") are precisely the sets V' x . [For example, we could connect the vertices 
of each V' x by a path or a cycle or a complete graph.] Now let G' = (V, E' U E"), and 
let Z(q,v) be the Potts-model partition function Zc(q, {v e }) with weights v e = — 1 
for e G E' and v e = v for e 6 It then follows immediately from (jl.2|) that 

\E"\ 

Z(q,v) = J2 vk P^) ( B - 9 ) 

fe=o 

with highest-order coefficient 

p m (q) = P G {q). (B.10) 

Now let G'(p) be obtained from G' by replacing each edge in E" by the 2-rooted 
graph 9^ s ' p ) (rooted at the endvertices). By construction, G'(p) is 2-degenerate, hence 
3-colorable. Moreover, G'(p) is connected if G is, and G'(p) is bipartite if G is and 
s is even. Finally, if G is embedded on a surface S, then G'(p) can also be taken to 
be embedded on S, by drawing the vertex (x, y) G V' x on the edge xy very near to 
x, and by choosing E" to be a path (or cycle) connecting the vertices of each V' x in 
cyclic order as one circles x. 

It follows immediately from Proposition 12.21 and equation ()2.32|) that 

P G >(p)(q) = Mq) PlE " l Z(q,v eff (q)) (B.ll) 

where v e jj(q) is given by (jB.3|) . In particular, for q G Q s we have v e g{q) — > oo as 
p — > oo, so that 

lim C s (g)- p l B "lp G , (p) (g) = P G (g) . (B.12) 

Moreover, the convergence is clearly uniform on compact subsets of Q s . 

Claim (i) then follows immediately from (a) using Hurwitz' theorem [SZ1 p. 262] 
and the fact that Pq is not identically zero. To prove claim (ii), choose e > small 
enough so that Pg(?o _ e ) an d Pc(qo + e ) have opposite signs, and choose po large 
enough so that C s (q ± e)~'^ E ^Pg'(p){Qo =t e ) has the same sign as Pg{% ± e) for all 
V ^ Po] then use the intermediate value theorem. ■ 

Remarks. 1. Thomassen |115l Theorem 3.9] considers the special case s = 2 and 
restricts attention to real q > 2. He establishes flB.12|) for integer q > 2 by counting 
colorings, and then asserts that "this also holds" for real q > 2. The assertion is 
correct but requires a bit of argument; Proposition 12.21 provides a convenient way of 
filling in the missing details, and determining the exact set Q s on which (jB.12|) holds. 

2. One consequence of this result is that the upper zero-free interval for 3-colorable 
(or 2-degenerate) planar graphs is the same as that for all planar graphs (unless the 
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largest chromatic roots have even multiplicity, which seems unlikely). In commenting 
on this corollary, Thomassen |115t p. 506] asserted that there exist planar graphs with 
real chromatic roots arbitrarily close to 4; but this assertion apparently arises from a 
misunderstanding of the Beraha-Kahane |H] theorem that 4 periodic x n free triangular lat- 
tices have complex chromatic roots arbitrarily close to 4. In fact I do not know of any 
planar graphs with real chromatic roots arbitrarily close to 4. A study of triangular- 
lattice strips [60] has thus far found chromatic roots up to ~ 3.64; and Baxter's work 
jS] suggests that sufficiently wide and long pieces of the triangular lattice will have 
real chromatic roots up to at least the Beraha number 5 14 = 4cos 2 (7r/14) « 3.801938 
(see also jHO] for further discussion). Finally, Douglas Woodall (private communica- 
tion) has found a 23-vertex plane triangulation with a chromatic root at ~ 3.811475. 
But I do not know of any planar graphs with chromatic roots larger than this. It 
is therefore conceivable that the upper zero-free interval for planar graphs extends 
below q = 4. I thank Carsten Thomassen and Douglas Woodall for correspondence 
concerning these questions. 
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